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Abstract. Fluctuation dynamics in magnetization plateaux are a relatively 
poorly explored area in frustrated magnetism. Here we use muon spin relaxation 
to determine the fluctuation timescale and associated field distribution width 
in the partial magnetization plateau of Ca3Co206. The muon spin relaxation 
rate has a simple and characteristic field dependence which we model and by 
fitting to the data at 15 K extract a fluctuation timescale r = 880(30) ps and 
a fleld distribution width A = 40.6(3) mT. Comparison with previous results on 
Ca3Co2 06 suggests that this fluctuation timescale can be associated with short- 
range, slowly fluctuating magnetic order. 
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1. Introduction 

When magnetic ions are arranged in certain structures their interactions can act to 
impede magnetic ordering. This is referred to as magnetic frustration and is widely 
studied [1] [5] . A well known example of a frustrated structure is antiferromagnetically 
coupled spins on a triangular lattice [3]. An alternative version of this is the 
antiferromagnetic coupling of ferromagnetic chains on a triangular lattice when the 
intrachain coupling is far stronger than the interchain coupling. This is realized in 
compounds such as CoNb2 06 t4j and Ca3Co206 where the moments within the 
chains have Ising (easy-axis) anisotropy. 

In applied magnetic fields frustrated magnets can undergo transitions between 
different magnetic structures. In between such transitions the magnetization of the 
sample stays fixed over a significant range of fields: a magnetization plateau. These 
show a simple fraction of the saturation magnetization such as 1/3 because the 
individual atomic spins are either aligned or anti-aligned with the external field. Such a 
partial magnetization plateau is observed in the intermediate, partially ordered phase 
of CaaCosOg between 10 and 25 K [i H [H [10]. A field of ~ 0.5 T is sufficient to 
enter the 1/3 magnetization plateau and then a further transition into the saturated 
magnetic state occurs at 3.6 T. Below '-^ 10 K a series of magnetization steps are seen 
in single crystals for fields applied along the c-axis, with pronounced hysteresis [TlllOj. 
That these magnetization steps occur in multiples of 1.2 T has led to the suggestion 
that they represent quantum tunnelling of the magnetization [SJ [7] 

Muon spin relaxation (/xSR) can be used to investigate both the static and 
dynamic properties of frustrated magnets [111 1121 113j . This is most commonly carried 
out in zero applied field, since the technique is rare in allowing this and is exquisitely 
sensitive to small magnetic fields that could emerge due to weak magnetic order. It 
is sensitive to fluctuating magnetic fields on time scales ranging from around 10~^^ 
to 10^^ s, intermediate between neutron scattering and ac magnetic susceptibility. 
Measurements on magnetic systems in applied fields along the initial direction of the 
muon spin allow the implanted muon spins to be decoupled from the field distribution 
they experience at their stopping site, providing a further window on the static and 
dynamic properties [12j . 

In this paper we describe the results of a /iSR investigation of the magnetization 
plateau in Ca3Co206. The broad and experimentally accessible plateau in the 
intermediate temperature region offers an ideal ground for quantitative comparison 
between models for the relaxation dynamics and experimental data. We propose a 
simple model to describe the field dependence of the muon spin relaxation rate, find 
that it describes the data successfully in the plateau region, and use it to estimate the 
time scale and field distribution associated with the magnetic fluctuations that the 
muon probes in this system. 

2. Experimental 

Stoichiometric Ca3Co206 has been synthesised by the conventional solid state reaction 
technique. High purity (> 99.99 %) starting materials CaCOs and C03O4 were dried 
at 150°C for 6 h before weighing. Mixed chemicals were calcined and sintered at 
810 — 860°C for 48 h with intermediate grinding. Finally the sintered powder was 
pressed into cylindrical rods and sintered at 925°C for 48 h in air. 

Our /j,SR experiment [Til [HI [13] was carried out on the newly constructed HiFi 
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Figure 1. Examples of the Pz{t) data recorded at different fields for Ca3Co206 
at 15 K. The fines plotted are fits of the data to Equation \3\ adding a slowly 
relaxing background at T. 

spectrometer [H] at the ISIS Pulsed Muon Facility. This has a 5 T superconducting 
magnet, almost ideally suited to the field range needed for this experiment. The 
powder sample was mounted in a 25 ^m silver foil packet on a silver backing plate 
and cooled to 15 K using a helium flow cryostat. Silver is used because the muon spin 
relaxation rate for muons stopping outside the sample will be close to zero over the 
field range of interest here. 

The fully spin polarized muons are implanted into the sample and generally stop 
at interstitial positions within the crystal structure without significant loss of spin 
polarization. The spin polarization Pzit) of the muon depends on the magnetic 
environment of the of the muon stopping site and can be measured using the 
asymmetric decay of the muon into a positron, which is detected, and two neutrinos. 
Around 10 million decay positrons were detected in each data set. The emitted 
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positrons are detected in scintillation counters around the sample position [T3] , divided 
into two banks forward (F) and backward (B) relative of the initial muon spin direction. 
The time-dependent asymmetry A{t) between the count rates in the banks of detectors 
is: 

^ ' NF{t) + aNB{t)' ^ ' 

where the parameter a describes the relative counting efficiency of the two banks of 
detectors, depending on the sample and detector geometry. From A{t) it is possible 
to infer the spin polarization Pz (t) of the implanted muon ensemble: 

P (t) - ~ ^'^g (2) 

The measured asymmetries A{0) and Abg are the initial and background values 
respectively and depend on the detector geometry and the fraction of the muon beam 
incident on the sample. The field variation of a, ^(0), and Abg required careful 
consideration in fields ^ 1 T because the decay positrons spiral significantly in the 
field. This was corrected for by comparison with data from reference samples |14| . 

Following on from earlier muon work on Ca3Co2 06 in the intermediate 
temperature region |8l , we expect that in close to zero applied field the muon spin 
relaxation will not show the full initial asymmetry observed in the paramagnetic region 
because the internal fields are too large to be observed at ISIS. (The pulsed nature 
of the muon beam gives a range of arrival times for niuons entering the sample which 
limits the maximum observed relaxation rate.) Takeshita and co-workers investigated 
the form of the spin relaxation up to around 1 T in this temperature region, finding 
that the initial asymmetry had been recovered by 0.5 T and that the relaxation was 
exponential over the whole field range. This agrees with the form of our data from 
up to 3.8 T, as shown in Fig. [T] and we therefore fitted our data to the function: 

=exp(-Ai), (3) 

where A is the relaxation rate. For A > 0.1 MHz, A can be extracted from the data 
independent from any field variation of a, ^(0), and A\^g. This is the case throughout 
our measured field range. In zero field A^g decreases slowly with time due to muons 
stopping in the silver sample holder and any weak fields perpendicular to the initial 
muon spin direction, but in applied field Abg is effectively constant. 

To model the field dependence of A we make the assumption that the distribution 
of fluctuating magnetic fields approximates to a Lorentzian [121 IH] with a width A 
and that one correlation time r is an adequate description of the fluctuation rate. We 
further assume that within the partial magnetization plateau these two parameters 
can be taken to be constant. Making these assumptions allows us to simplify Rcdficld's 
equation [15] to the form: 



The parameter 7^ = 27r x 135.5 MHzT"-'^ is the muon gyromagnetic ratio 



3. Results 



In Fig. [5] (a) we show the relaxation rate A extracted from the data at 15 K as a 
function of the magnetic field Slf = t^oH. Three behaviours are apparent in the 
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Figure 2. Parameters extracted from the /iSH data using Equation [S] (a) 
Relaxation rate A plotted against the applied field hqH. (b) Inverse of the 
relaxation rate plotted against the square of the applied field /Xgi?^. The 
trend obtained by fitting the data to Equation |4] is shown in the shaded 
magnetization plateau region. 



data: At fields below 0.5 T there is a rapid decrease in A, between 0.5 and 3.6 T 
the decrease in A with increasing field is slower, and above 3.6 T any change in A is 
smaller than the experimental error. These regions correspond to the initial increase 
in bulk magnetization, the magnetization plateau, and the saturated magnetic phase 
respectively. 

To make a more obvious comparison with the Redfield equation Q discussed 
above we plot A~^ against ^i^H^ in Fig. [2] (b). This rearrangement means that 
the trend for constant A and r will appear as a straight line, as is observed in the 
data in the magnetization plateau region highlighted. Fitting equation 0] to the data 
shown in Fig. [5] between 0.5 and 3.6 T gives the parameters: A = 40.6(3) mT and 
T — 880(30) ps. This gives the line plotted in Fig. [2] (b). These parameters are 
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insensitive to small reductions in the fitting range at either end of the plateau region. 

4. Conclusion 

Using the partial magnetization plateau in Ca3Co206 as a model system in which 
to test a simple model for the field-dependence of the muon spin relaxation rate in 
a magnetization plateau we find very good agreement between the form predicted 
by the model and the experimental data. We are able to obtain a field distribution 
width and fluctuation timescale for the dynamic magnetism within the plateau. This 
suggests that this model could also be applied to other magnetic systems that exhibit 
a magnetization plateau and provide new information about their dynamic behaviour. 

The width of the field distribution that we obtain is comparable with the 
~ 100 mT static internal field observed in zero- field measurements [5|. However, 
were the muon probing such static fields they would have been decoupled from the 
muon spin at fields lower than those in the magnetization plateau. Both the previous 
/xSR results [81 and neutron diffraction work [10 show strong evidence for shorter 
ranged, dynamic correlations coexisting with the long-range magnetic order at this 
temperature. Such dynamic correlations would not be decoupled from the muon spin 
so easily and it is therefore this behaviour that we are probing in this experiment. 
Since the field distribution width will be determined by the same magnetic moments 
the same distance from the muon stopping site(s), and this is entirely consistent with 
our results. 

The timescale we find for the fluctuations < 1 ns sits comfortably within the 
sensitivity range of /iSR and suggests that these fluctuations would be static on the 
timescale of neutron diffraction. Again, this is consistent with the previous neutron 
diffraction measurements where the short-range magnetic order is quasistatic [TU] , and 
the previous /iSR results where a significant strongly-relaxing component suggestive of 
slowing fiuctuating short-ranged order is observed at this temperature [8l[9j. Earlier ac 
susceptibility measurements found a significantly longer timescale, ~ 0.4 s [7], which 
is far outside the window we can probe using muons, and is likely to be associated 
with a separate process within the system. 

Comparing our model to the data in the region B < 0.5 T we find far poorer 
agreement. Of course, the assumption that A is independent of field is unlikely to 
be appropriate here and we would only get some approximation to the average value. 
The parameters coming from our data for B < 0.5 T are: A ~ 30 mT and t ~ 5 ns, 
which suggest that the decoupling from longer ranged order remains relevant in this 
field range. Takeshita et al. carried out a similar analysis of their data at 20 K and 
obtained larger, although comparable, values for both parameters [9]. 
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